Detecting correlation functions of ultracold atoms through Fourier sampling of 

time-of-flight images 
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We propose a detection method for ultracold atoms which allows reconstruction of the full one- 
particle and two-particle correlation functions from the measurements. The method is based on 
Fourier sampling of the time-of-flight images through two consecutive impulsive Raman pulses. For 
applications of this method, we discuss a few examples, including detection of phase separation 
between superfluid and Mott insulators, various types of spin or superfluid orders, entanglement, 
exotic or fluctuating orders. 
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Ultracold atomic gas provides an ideal system for 
controlled study of various kinds of strongly correlated 
many-body physics Q, |^ . To reveal different many- 
body phenomena, the correlation function plays a criti- 
cal role. In condensed matter systems, information about 
the correlation function is typically collected through lin- 
ear response or scattering experiments 0. For ultracold 
atomic gas, the most powerful detection technique is ar- 
guably the time-of-flight imaging For ballistic ex- 
pansion of the atomic gas, the time-of-flight images ac- 
tually give information of the atom number distribution 
in the momentum space, which corresponds to the diago- 
nal terms of the one-particle correlation function. As we 
have only diagonal correlation function in the momentum 
space, it is in general inadequate to use it to reconstruct 
the correlation function in the real space. 

In this paper, we propose a detection method which 
extends the time-of-flight imaging technique by gaining 
much more information of the correlation function. The 
method detects all the non-diagonal correlation terms as 
well as the diagonal ones in the momentum space through 
application of two consecutive impulsive Raman pulses 
at the beginning of the expansion which introduces a 
tunable momentum difference to the correlation terms. 
With this kind of Fourier sampling in the momentum 
space, one can reconstruct the full one-particle correla- 
tion function in the real space. If one repeats the time- 
of-flight imaging measurements many times, as has been 
shown in recent works 0, it is also possible to de- 
tect the noise spectroscopy by looking at the statistical 
correlation between different images. With a combina- 
tion of the Fourier sampling and the noise spectroscopy, 
one can also reconstruct the full two-particle correlation 
function. The reconstruction of the full one-particle and 
two-particle correlation functions represents an unprece- 
dented detection ability, which gives very detailed infor- 
mation of the underlying many-body system. As illustra- 
tion of applications, we discuss a few examples, including 
detection of phase separation between the superfluid and 
the Mott insulator states for bosonic atoms in a trap, var- 
ious types of spin or superfluid orders in multi-component 
Bose or fermi gas, entanglement in an optical lattice, pat- 
terns of the valence bonds associated with some exotic 



quantum phases , and the fluctuating orders which ex- 
ist only in short distance p^ . 

In the time-of-flight imaging experiment, what one 
measures from light absorption is the column integrated 
density of the expanding atomic cloud. The signal is 
proportional to the column average of the atom density 
operator, which, for the spin component a, is denoted 
as (n„ (r,i)) = (*t (^r,t) (r,t)), where (r,t) repre- 
sents the field operator for bosonic or fermionic atoms 
at the position r and the expansion time t. We assume 
ballistic expansions with the atomic collision effect negli- 
gible during the time of flight HtI. In such a case, if the 
expansion time t is long enough so that the size of the 
final expanded cloud is much larger than the size of the 
initial one, the density (tIq, (r,t)) is connected with the 
initial momentum distribution of the atoms by a simple 
relation (n^ {r,t)) cx {ua (k)) = (vJ/J^ (k) vj/^ (k)) with the 
corresponding wave vector k —mr/ {fit) ^J. So the 
conventional time-of-fiight imaging measures the diago- 
nal terms of the one-particle correlation function in the 
momentum space. 

To reconstruct the full one-particle correlation func- 
tion, it is required to measure also the non-diagonal cor- 
relation terms (4' J, (k) ^q, (k')) in the momentum space. 
For that purpose, we propose a detection method as il- 
lustrated in Fig. 1, which combines the time-of-flight 
imaging with two consecutive impulsive Raman pulses. 
We assume there is an additional atomic hyperfine spin 
level [3 which is initially empty (such a level is always 
available for typical alkali atoms). Right after turnoff 
of the trapping potential but before any significant ex- 
pansion of the atomic cloud, we apply two impulsive Ra- 
man pulses to all the atoms. The duration 5t of the 
Raman pulses is short so that one can neglect the cloud 
expansion and the atomic collision within St. We as- 
sume that the two travelling-wave beams for the first 
Raman operation are propagating along different direc- 
tions with the corresponding wave vectors ki and k2, so 
the effective Raman Rabi frequency has a spatially vary- 
ing phase with the form O (r) = r2oe'*^'''' ''+'^i\ where 
(5k = k2 — ki and (pi is a constant phase. The Hamil- 
tonian within the interval 5t can then be written as 
H = / d^rl^ (r) ^E'Jj (r) (r) + H.c. Transferring this 
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Hamiltonian into the momentum space, we have 
H = J2 f^oe^'^^^i (k) (k+,5k) + H.c. 

k 



(1) 



We choose the intensity of the Raman beams so that 
il.f)dt = 7r/4. After this Raman operation, the final field 
operators in the momentum space, denoted as ^'^ (k) and 
(k) , are connected with the initial ones (k) , vjf^ (k) 
through the relation 

<(k) = [*„(k) + e*'^i*^(k+5k)] /n/2, 

^'^{k) = [v|/^(k)-e-'^ivl/„(k-5k)]/\/2. (2) 

Then, we immediately apply the second Raman opera- 
tion which is from two co-propagating laser beams con- 
necting the levels a, (3 with a spatially constant effec- 
tive Rabi frequency fioe*'''^- With the same pulse area 
n^St = 7r/4, the second Raman operation induces the 
transformation 



<(k) = [^'^{k) + e^^^%{k)]/V2, 

^'^{k) = [%{k)-e-^^^K{k)]/V2, (3) 

where 5*'^ (k) , 'i>'^ (k) denote the final field operators af- 
ter the second Raman operation. We then perform 
the conventional time-of-flight imaging, with the atoms 
in different spin components a,/3 separated during the 
flight through a magnetic field gradient 5]. This imag- 
ing measures the momentum distribution {n'^ (k)) = 

(4"^t(k)vl/'^(k)) and (n'^ik)) = (*^t (k) M/^ (k)). 

With a combination of Eqs. (2) and (3), we find that 

the difference between the two images ( n'^^ (k) \ = 



(n^ (k)) — yi'^ (k)y gives the non-diagonal correlation 
terms 

«^ (k)) = - Re (k) {k-5k) e^*>) , (4) 

where dip = ip2~'^i a-nd we have used the fact that ^'/j (k) 
is initially in the vacuum state. 
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FIG. 1: (lA) Raman pulses with different propagating di- 
rections introduce a momentum kick. (fB) The momentum 
ladder connected by two consecutive impulsive Raman pulses. 
The first 7r/2 Raman pulse (pulse f ) introduces a tunable mo- 
mentum kick 5k, while the second one (pulse 2) has no kick. 
The atoms are initially in the spin-Q state. 

With the above method, we have shown how to 
measure the real and imaginary parts of the non- 
diagonal correlation function (^"^^ (k) ^'^ (k— (5k)) by 



choosing the relative phase Sip = or 7r/2, respec- 
tively. The momentum difference 5k is controlled 
by the relative angle of the two laser pulses for the 
first Raman operation. By varying this angle, 5k 
can be varied from to 2ko, where |ko| — 27r/A is 
the optical wave vector. With a Fourier transform of 
(4' J, (k) (k— (5k)), one can reconstruct the real-space 
correlation function (vpt (r+(5r/2) {^-5v/2)) = 
J {^l (k+(5k/2) (k-,5k/2)) e*(''-*'-+'5k--)d5kdk/27r, 
with a spatial resolution in r down to 27r/ l^kj^^^^^ = A/2. 
This kind of Fourier sampling in the momentum space 
through the laser phase gradient allows us to directly 
probe very small spatial structure in the ultracold atomic 
gas, although the probe beams are always shincd on all 
the atoms without separate addressing of any particular 
region. For instance, for the case of atoms in an optical 
lattice, with such a resolution, we can reconstruct the 
correlation (aj, (i) Qa (j)) of the mode operators at 
two arbitrary sites i and j (the lattice spacing is A/2). 

Before showing how to measure the full two-particle 
correlation function, we add a few remarks here. First, 
note that in the above method, it is essential to apply two 
consecutive Raman operations. If we introduce momen- 
tum difference 5k by applying only one Raman operation 
on the same spin component a with the Hamiltonian 
H' = J2k^oe''^'-i'i{k)^aik+6k) + H.c, the H' will 
couple the whole momentum ladder '^a (k), '^a {k±5k), 
(k ± 25k), • • •, and we thus cannot get a simple form 
of transformation such as those given by Eqs. (2) and 
(3). So the method with two consecutive Raman opera- 
tions seems to be the easiest one for measuring the non- 
diagonal momentum correlations. Second, in the above 
method, we have assumed to tune the momentum differ- 
ence by changing the relative angle of two laser pulses. 
This tuning can also be achieved by applying a sequence 
of laser pulses incident along two fixed directions with 
a small angle 59, with each pair of pulses introducing a 
fixed momentum kick k^ = 2 |ko| sin {56/2) ~ |ko| 59 (see 
illustration in Fig. 2 and its caption) . With a maximum 
of N such pulse pairs, the momentum difference 5k in 
the correlation function (4) can be tuned among different 
values k^, 2ks, • • •, A^k^, which corresponds to a discrete 
Fourier sampling of the atomic correlation function with 
the momentum-space and real-space resolutions given re- 
spectively by I kg I (56' and X/ {N5d). Finally, the above 
method can be generalized straightforwardly to measure 
also the spin-spatial correlations. If we have multiple Zee- 
man spin components a,a' , ■ ■ ■ for the state of the atomic 
cloud, one can reconstruct the full spin-spatial correla- 
tion (*t (ki)*„, (ka)) (or (*t (r,) (rs))) through 
a combination of the above Fourier sampling with a pair 
of Raman pulses which mixes the spin components a, a'. 

The above method, combined with the 
noise spectroscopy, can also be used to recon- 
struct the full two-particle correlation function 
(*a(i"i)'^(r2)*a(r3)*a(r4)). For the noise spec- 
troscopy la, 13 1 one just needs to note that for different 
runs of experiment, the time-of-fiight images can have 
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FIG. 2: Tunable momentum kicks with pulses propagating in 
fixed directions. Steps (A) and (D) represent 7r/2 pulses with 
no momentum kicks. Steps (B) and (C) represent tt pulses 
which a small fixed momentum kick ks (pulses for B, C are 
propagating in reverse directions). The intermediate pulses 
B, C can be applied with a controlled number of times, which 
correspond a discrete Fourier sampling of the correlation func- 
tion with a momentum-space resolution of ks- 



quantum fluctuation even if one starts with the same 
state of the atomic cloud. One image corresponds to a 
single-run measurement of the column integrated atomic 
density operator, and by looking at correlation of differ- 
ent images, one can find out quantum correlation of the 
density operator. To detect such quantum correlation, 
all the technical noise for the images needs to be reduced 
below the level of quantum noise. It is remarkable that 
the recent experiments have shown that such quantum 
correlation of images is indeed visible after suppression of 
the technical noise 7] . Each time-of- flight imaging corre- 
sponds to a measurement of the one-particle momentum 
correlation (^J^, (k) vJ^q, (k — (5k)) . To reconstruct the 
full two-particle correlation function, one just needs to 
repeat each of such imaging measurements M times (M 
needs to be sufficiently large so that 1/\/M, which char- 
acterize the statistical error for the quantum correlation 
measurement, is sufficiently small). Then, by looking at 
the correlation of arbitrary two images corresponding 
to (vft (k)vl/„(k-5k)) and (vf^ (k') (k' - (5k')) 
respectively, one gets the two-particle correlation 
(4'l(k)^'„(k-,5k)*t^(k')«'a(k'-5k')). As the 
two images can be arbitrarily chosen (with different 
k, k', (5k,(5k'), this gives us the full two-particle cor- 
relation in the momentum space, and the real-pace 
correlation (^-t (rj) (ra) ^-J, (rg) (r4)) is simply 
its Fourier transform. The spatial resolution in this 
case is still given by A/2 similar to the one-particle 
correlation measurement. 

The measurement of the full one-particle and two- 
particle correlation functions gives very detailed infor- 
mation of the many-body properties of the underlying 
system. Here, we give a few examples for illustration 
of its applications. We start with one-particle correla- 
tion function. As the first example, note that for the 
superfluid to Mott insulator transition in a global har- 
monic trap (as it is the case for experiments), theory 
has predicted phase separation with layers of superfluid 
states intervened in by layers of Mott insulator states 
with decreasing integer filling numbers Q . Such a picture 



have not been confirmed yet by experiments as with the 
conventional time-of-flight imaging, it is hard to see this 
layer-by-layer structure. Through the Fourier sampling, 
however, one can reconstruct the one-particle correlation 
(^'jj (ri) "ifa (r2))- From this correlation, it would be evi- 
dent to see the phase separation: the Mott insulator state 
is characterized by a vanishing non-diagonal correlation 
and a constant (integer) diagonal correlation, while the 
superfluid state is characterized by both non-zero diag- 
onal and non-diagonal correlations which vary continu- 
ously in space. 

With ultracold atoms in an optical lattice, one can re- 
alize different magnetic Hamiltonians 12]. Depending on 
the lattice geometry and interaction conflgurations, such 
Hamiltonians may support various types of magnetic or- 
ders 0, . The magnetic orders typically can be writ- 
ten in the form (Si) = ViCOs(K-ri) -f V2 sin (K • r^) 
0, where is the spin operator on the site i with 
the coordinate r^, Vi and V2 are two vectors specify- 
ing certain directions, and K characterizes spatial vari- 
ation of the order parameter. The spin Si is deflned 

as Si — J2aa' flLfaa'flia'/2, whcrc CT is the Pauli matrix 
and aia is the mode operator which is connected with the 
fleld operator 'i'a (r) through aia = / w* (r) (r) (Pr 
(w* (r) is the Wannier function for the site i). Through 
Fourier sampling of the time-of-flight images with a res- 
olution A/2, one can reconstruct correlation of the mode 

operators (al^aja'^ on any sites {i — j for the special 

case) . So one can directly detect any type of spatial vari- 
ation of the order parameter (Si) (note that the con- 
ventional time-of-flight imaging cannot detect spatially 
varying order parameters). 

Another interesting application of this Fourier sam- 
pling technique is to detect local fluctuating orders. Lo- 
cal fluctuating orders, such as stripes (unidirectional den- 
sity waves), have received wide attention in strongly cor- 
relatedphysics (in particular for high-Tc superconduc- 
tors) [T3|. The local fluctuating order typically takes 
place near a critical point with competing orders or at the 
proximity of an ordered phase. With ultracold fermions 
in an optical lattice, the fundamental Hamiltonian is 
very similar to those models for high-Tc superconduc- 
tors [13 . Hsj . One expects that competing fluctuating 
orders may arise as well in the phase diagram of this 
system. The local fluctuating orders such as stripes cor- 
respond to real-space patterns of some micro-phase sep- 
aration p^. and the conventional time-of-flight imaging 
cannot see it because it will be averaged out. However, 
through the Fourier sampling technique with a spatial 
resolution down to the lattice constant, it should be ev- 
ident to detect these local fluctuating orders whenever 
they show up. A measurement of the one-particle cor- 
relation ('I'Jj (ri)^'^/ (r2)) is typically enough to probe 
the fluctuating orders such as the local density waves 
(stripes). 

The above discussion shows some applications of the 
one-particle correlation. With the two-particle correla- 
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tion measured, one can gain further information. For ex- 
ample, one can use the two-particle correlation to unam- 
biguously detect the entanglement pattern for ultracold 
atoms in an optical lattice. Through controlled atomic 
collisions, one can generate entanglement for atoms on 
different lattice sites and some initial experimental 
evidence of entanglement has been reported for such a 
system It is hard however to directly detect atomic 
entanglement between different sites or regions inside an 
optical lattice as the conventional detection technique 
donot have such a spatial resolution. This entanglement 
can be unambiguously confirmed with a measurement of 
spin correlations for atoms on different sites (similar to 
the Bell inequality measurement). As each spin operator 
involves only two atomic mode operators, the spin cor- 
relation is included in the two-particle correlation func- 
tion. So, a measurement of the two-particle correlation 
can gives information of spin entanglement pattern in the 
lattice. 

For fermionic systems with pairing instability (such 
as the fermionic superfluid state 0), it is also desirable 
to have a method to directly measure the Cooper 
pair function (^'^ (ri) (r2)) . The two-particle cor- 
relation also gives information of such pair function 
as with Cooper pairing, the two-particle correlation 
can be typically approximated with the following 
decomposition (^^l (ri) ^f^, (rs) (rs) ^' 

(*t (n) (r2)) (VI/,. (rg) (r4)) + (ri) (r4)) 

{K' (1-2) ir3))-{H (ri) (rs)) {^l (r^) (r4)) 
by applying the Wick theorem. As the one-particle 
correlation is known already with the Fourier sampling, 
so a measurement of the two-particle correlation gives 
the pair function {"if a {^i)'^ a' {^2)) (the two-particle 
correlation actually contains more information, for in- 
stance, it can be used to check self-consistently whether 
the above decomposition is valid). The pair function 



{"^a (ri)^'Q' {1C2)) can tell us the symmetry (s-wave or 
d-wave for instance) and the size of the Cooper pairs as 
well as how the pair structure changes in space. 

Finally, a potentially more interesting application of 
the two-particle correlation might be that it gives us a 
way to directly detect the patterns of valence bonds in 
an optical lattice. A valence bond on the sites i,j (not 
necessarily neighbors) is simply the ground state of the 
bond operator Qij = Sj • Sj (Si are spin operators) 0. 
The resonating valence bond (RVB) states P,'!^, which 
include many possible patterns of the valence bond dis- 
tribution in the lattice, have been conjectured as one of 
the most likely ground state of the t-J or Hubbard mod- 
els in the strongly correlated limit. As the number of 
possible configurations of the RVB states increses expo- 
nentially with the size of the lattice, it is hard to figure 
out the distribution pattern of the valence bonds in an 
lattice. However, with an atomic realization of the t-J 
model p^ . one can directly measure the bond opera- 
tors Qij in an optical lattice to find out the most likely 
distribution pattern of the valence bonds. Each bond 
operator Qij corresponds to a special component of the 
two-particle correlation function, so a measurement of 
the two-particle correlation with Fourier sampling gives 
complete information about the bond distribution. 

In summary, we have proposed a detection method 
which combines the commonly-used time-of-flight imag- 
ing technique with the Fourier sampling based on appli- 
cation of two consecutive impulsive Raman pulses. This 
detection method allows us to reconstruct the full one- 
particle and two-particle correlation functions. We have 
discussed a few examples for illustration of the wide ap- 
plications of the correlation functions for probing many- 
body properties of the underlying system. 

This work was supported by the NSF award (0431476), 
the ARDA under ARO contracts, and the A. P. Sloan 
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